Occupation of a resonant level coupled to a chiral Luttinger liquid 
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We consider a resonant level coupled to a chiral Luttinger liquid which can be realized, e.g., at a 
fractional quantum Hall edge. We study the dependence of the occupation probability n of the level 
on its energy e for various values of the Luttinger-liquid parameter g. At g < 1/2 a weakly coupled 
level shows a sharp jump in n(e) at the Fermi level. As the coupling is increased, the magnitude 
of the jump decreases until \/2g, and then the discontinuity in n(e) disappears. We show that n(e) 
can be expressed in terms of the magnetization of a Kondo impurity as a function of magnetic field. 

PACS numbers: 71.10.Pm, 73.43.Jn, 73.21.Hb, 72.10.Fk 
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The transport of electrons in low-dimensional struc- 
tures is strongly affected by electron-electron inter- 
actions. This effect is particularly strong in one- 
dimensional (ID) systems, where the interacting elec- 
trons form the so-called Tomonaga-Luttinger liquid jj]. 
The main feature of this system is the power-law sup- 
pression of the tunneling density of states at low energies. 
This effect has been recently observed experimentally in 
ID electron systems realized in quantum Hall edges [||, 
carbon nanotubes || [| , and quantum wires 0] . 

A generic property of mesoscopic devices is the pres- 
ence of imperfections that affect the transport of elec- 
trons. The most important are the defects that form 
energy levels near the Fermi level and give rise to strong 
scattering of electrons. In this paper we consider a single 
level coupled to a ID system exhibiting Luttinger liq- 
uid behavior. Experimentally the level can be formed 
either by a random impurity or by an artificially created 
quantum dot in the vicinity of the ID system. In recent 
experiments J5|, |(| resonant tunneling via such a level 
has been observed. The theory of resonant tunneling via 
a level coupled to two ID conductors was developed in 
Refs. § || . 

In contrast to the previous theoretical work |?], §L |), 
we study a resonant level coupled to a single ID con- 
ductor. The main quantity of interest is the occupation 
probability n of the level as a function of its energy e 
(measured from the Fermi level). Clearly, at large pos- 
itive energy s the level is empty, n = 0, and at large 
negative energy it is filled, n — 1. If the level is cou- 
pled to a Fermi-liquid lead, the occupation probability 
changes continuously from to 1 when the energy e is 
within the level width V from the Fermi level: 



Hamiltonian H = Hq + H t , where 
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We show below that the dependence n(e) is modified sig- 
nificantly in the case of coupling to a Luttinger-liquid 
lead. Most importantly, if the interactions are strong 
enough, n(e) has a discontinuity at e = 0. 

We model the level coupled to a ID conductor by the 
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The first term in Eq. (g) describes the resonant level 
with energy e; the fermion operator creates an elec- 
tron in this state. The second term in Eq. ^ describes 
the Tomonaga-Luttinger liquid in the lead in terms of 
a chiral boson field f(x) with commutation relations 
[ip(x), f(y)] = i7rsgn(x — y). The chiral Luttinger liquid 
model is the most natural description of fractional quan- 
tum Hall edge states pfj[ . In the case of a non-chiral ID 
system coupled to the level at only a single point x = 
one can neglect the "odd" bosonic modes and convert 
Jll| the Hamiltonian to the chiral form (||). The Hamil- 
tonian Ht describes the tunneling of the electron between 
the level and the point x = in the ID system. We use 
the standard bosonized expression for the annihilation 
operator ip of an electron in the chiral Luttinger liquid in 
terms of the field f(x), the Luttinger-liquid parameter g, 
the bandwidth D and the Fermi velocity v. In this paper 
we only consider the case of spinless electrons. Many of 
our results, including the existence of the discontinuity 
in n(e) at g < 1/2, are not sensitive to this assumption. 
Also, in the case of tunneling into edge states of quan- 
tum Hall system with filling factor g < 1 all electrons are 
spin-polarized by the magnetic field. In other cases the 
level can be polarized by applying the magnetic field. 

We first calculate the occupation probability n(e) in 
the second-order perturbation theory in H t . At e > 
the result has the natural form 
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Parameter T = t 2 /2hv has the physical meaning of the 
level width in the Fermi- liquid case g — 1. The tunneling 
density of states v g (bS) of the Tomonaga-Luttinger liquid 
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has the well-known form 
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where T(x) is the Gamma function. Due to the electron- 
hole symmetry of the problem, the occupation probabil- 
ity at negative energies can be found as n(e) = 1 — n(— e). 
At g = 1 the expression (Q) coincides with expansion of 
the Fermi-liquid result (|l|) at e 3> I\ 

At g > 1/2 the integral in Eq. (Q) diverges at e — > 0. 
This means that in order to find n(e) near the Fermi level 
one has to sum up the whole series of the perturbation 
theory in t. On the other hand, at g < 1/2 the integral 
is convergent at all e, i.e., for a weakly coupled level the 
second-order perturbation theory is sufficient to describe 
n(e) at any energy. For instance, at e — ► +0, we find 
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At e <C D the occupation probability approaches its limit 
(0) following a power-law dependence: 
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Since n(+0) ^ n(— 0) = 1 — n(+Q), the system reaches 
two different ground states at e — > ±0, and the occupa- 
tion probability n(e) experiences a finite discontinuity at 
the Fermi level: A = n(— 0) — n(+0). The divergence 
of the perturbation theory at g > i indicates that in 
this regime the coupling of the level to the lead is a rele- 
vant perturbation that should lift the degeneracy of the 
ground state at e — 0, resulting in continuous n(e). This 
can be shown formally using the mapping to the Kondo 
problem discussed below. 

As the coupling T of the level to the ID system grows, 
the magnitude of the jump A = 1 — 2n(+0) at g < | 
decreases. We will now show that at a sufficiently strong 
coupling the jump disappears. Let us consider a semi- 
infinite quantum wire in the region x > with a scat- 
terer at point x = L. The system is described by the 
Tomonaga-Luttinger model with a scatterer: 
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dx — w cos(2ttN). 

(8) 

Here N — [kpL + 9(L)]/ir is the number operator of 
electrons in the region < x < L. The fields 8 and <j> 
satisfy the commutation relation [9(x), 4>(y)] = (z7r/2)[l + 
sgn(x — y)] and the boundary condition Q(x = 0) = 0. 

If the scatterer is very strong, w — ► oo, it effectively 
cuts the region < x < L from the wire. This region be- 
comes a quantum dot, coupled to the rest of the wire by 




FIG. 1: Schematic phase diagram for the behavior of the oc- 
cupation probability n(e). The discontinuity occurs at g < | 
and weak enough effective coupling F/D. 



tunneling through the barrier created by the scatterer. 
At large w the field 9 at point x = L is pinned at one 
of the values 0(L) = nm — kpL, where integer m is the 
number of electrons in the dot. The ground state energy 
of the dot is then E m = (hv /2-KgL)(nm — kpL) 2 . For 
a given L the number of electrons in the dot that mini- 
mizes E m is the integer nearest to kpL/n. One can easily 
check that at kpL/ir = M + I — -j^j£ with integer M 
and positive s the ground state of the dot has m = M 
electrons, and the energy of adding (M + l)-st electron 
is Em+i — Em = e. Thus at large w the Hamiltonian (|^) 
is equivalent to the model (|) and of a resonant level 
coupled to a ID system pi . 

One can now increase the coupling by lowering w. 
In the limit of very strong coupling, w — > 0, the scat- 
tering term in the Hamiltonian (|8|) can be accounted 
for in first-order perturbation theory. In an infinite 
ID system the thermal average of the scattering poten- 
tial shows power-law behavior at low temperatures as 
(T/D) 9 w cos(2kpL). However, if the distance L from the 
scatterer to the end of the wire is finite, the renormaliza- 
tion of the scattering potential is cut off at the low energy 
scale equal to the level spacing in the dot, ~ hv/L. Al- 
though the renormalization of the scattering potential in 
a semi-infinite wire may be large, it remains finite, and 
at sufficiently small w one can use the first-order result 
5Eq ~ (hv/LD) 9 wcos(2kpL). Thus at strong coupling 
the energy is an analytic function oikpL which shows no 
singularities at degeneracy points kpL = ir(M + ^) where 
e = 0. The analyticity of E (e) means that the ground 
state is non-degenerate, and the occupation probability 
n(e) = dE /de (see Eqs. (||) and (||)) is continuous. 

From the perturbative calculations in the weak- and 
strong-coupling limits we deduce a schematic phase dia- 
gram, Fig. [I]. The occupation probability n(e) is a con- 
tinuous function of e at g > h (including the Fermi- liquid 
case g = 1) and also at g < g , if the coupling to the lead 
is sufficiently strong. On the other hand, if g < h and 
the coupling is weak, n(e) has a discontinuity at e = 0. 
As the coupling is increased at g < i, the discontinuity 
A decreases until the phase boundary is reached, above 
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which A = 0. We now show that A is not a continuous 
function of the coupling strength, and that just below 
the boundary A = ^J2~g. 

Let us consider a level relatively weakly coupled to the 
lead with g < i. The presence of a jump in the occu- 
pation probability n(e) means that at e — the ground 
state of the system is degenerate. One could expect that 
tunneling term (||) that couples the two degenerate states 
by moving an electron from the level to the lead would 
lift the degeneracy. This does not always happen because 
the tunneling of an electron into a Tomonaga-Luttinger 
liquid is suppressed at low energies, leading to the renor- 
malized level width T(e) oc e 7 . Thus at 7 > 1 the cou- 
pling is irrelevant, as it is much smaller than the level 
separation e. In the opposite case of 7 < 1 the width is 
much larger than the level separation at e — ► 0, so the de- 
genaracy of the ground state is lifted. Thus at the phase 
boundary in Fig. |l] the exponent 7=1. 

At t — ► the suppression of tunneling manifests it- 
self in the power-law behavior of the density of states 
(|J). Thus the exponent 7 = - — 1, and the condition 
7 > 1 of the degenerate ground state reproduces our 
earlier result g < |. At t 7^ the virtual hopping pro- 
cesses smear the quantization of the charge, and at low 
energies the perturbation Ht transfers only a fraction 
A = n(— 0) — n(+0) of an electron charge into the lead. 
Mathematically this can be accounted for by a renor- 
malized coupling term Ht <x t[a) 'e J¥, (°) A/ + h.c], c.f. 
Eq. (||). Thus the exponent of the tunneling density of 
states becomes 7 — ^ — 1 when only a fraction A of the 
electron charge is transferred. At the phase boundary 
the condition 7 = 1 is then equivalent to A = \f2g~. 

The phase diagram Fig. [j] and the critical value of the 
jump A = \/2g at the phase boundary are the main 
results of this paper. 

One can also study how 71(e) approaches the discon- 
tinuity at e = anywhere below the phase boundary in 
Fig. [j]. The discussion is essentially identical to the one 
leading to Eq. (R) at t — > 0. By substituting v{ui) oc u 1 
instead of (|1) into (|]), we find at small positive e, 

- n(e) OC e n 1 in{(A 2 /9)-2,l} for A > ^ (g) 

Additional support for the above results is found by 
mapping our model to two well-known problems: a two- 
level system with Ohmic dissipation and the Kondo 
problem. Let us define a unitary operator U\ — 
exp[iip(0)S z j \/g}, where S z — aJa — 1/2. Unitary trans- 
formation of the Hamiltonian (||) and (^) yields the 
Hamiltonian of the two-level system 

UIHU, =H + tJ^Ls* -*?S.*m + const . (10) 
V nnv y/g ax 

where S x — (a + er)/2. The second term describes hop- 
ping of a particle between the two levels, S z = ±i, and 



the third term is the coupling to harmonic oscillators with 
dissipation strength a=\/2g |l2). It is known that there 
is a quantum phase transition between a phase where a 
particle is localized in one of the two levels and another 
phase where it is delocalized. In the limit i — > the 
transition occurs at a = 1 Jl2] |, in agreement with our 
weak-coupling result Our analysis suggests that as 
the coupling t is increased in the localized phase, the 
difference A of the averages (S z ) over the two ground 
states at e = will decrease from 1 to ^/2g = a^ 1 / 2 , 
after which the system will suddenly delocalize, A = 0. 
To the best of our knowledge, only a brief discussion |TJ] 
can be found in the literature on the e dependence of the 
ground-state average (S z ) at a — 1 — 77, and the result is 
consistent with A = oT x l 2 to linear order in ?/. 

To understand the connection to the Kondo prob- 
lem, let us introduce another operator U2 = 
U\ exp[— i^/2ip(0)S z }. Unitary transformation yields 

UtHUi = H + tj-^r- ( S + e*^(°) + S-e-^vW) 
V 2-kTiv \ I 

* V2"VM2)+ ^nst, (11) 

where S + = and S~ = a. Equation ( |Tl| ) is also ob- 
tained by bosonizing the anisotropic Kondo model de- 
scribing electrons scattered by a spin-^ impurity, 

H K = J e k a} k ^a kll dk + J J J^a^al^^dkdk' + eS z , 

where e k — hvk, a 1 is the Pauli matrix, and summa- 
tion over repeated indices is assumed (/j, v =T,| and 
i = x,y, z). The Kondo couplings Ji are given by 

j,= t ^, Jz = 2^(1 --y, ( i2) 

where J x = J y = J]_. The occupation probability n{e) is 
directly related to the magnetization M. of the impurity 
spin S coupled to the Zeeman field e: n(e) = A4(e) + l/2. 
The phase diagram Fig. ^ is easily understood from the 
renormalization-group flow diagram of the Kondo prob- 
lem. The phase of a smooth n(e) curve corresponds to 
the antifcrromagnetic Kondo coupling, J z > — |Jj_|- In 
this case the impurity spin is completely screened as the 
magnetic field e is turned off: A4(±0) = 0. At g — h <C 1 
in the weak-coupling limit, T <C D, one can easily obtain 
n(e) from the Bcthe ansatz results [[[5] for A4(e). 

The phase of a discontinuous n(e) curve in Fig. |l| cor- 
responds to the ferromagnetic Kondo case, J z < —\J±_\, 
where the spin flip coupling Ji is rcnormalized to zero. 
One might then naively expect full moment to appear: 
A4 = — sgn(e)/2. This is not correct, however, as we have 
seen in Eqs. (g) and (0). In the weak-coupling regime we 
may use one- loop scaling equations for ji — Ji/(2Trhv), 
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where dl = — dh\D. For < ~j z we integrate them to 
I = ln(D/e) with initial condition M. = —1/2. We find 
2M(e) = - exp(j z - j*), where 



it _ J z -c+(c + Jz )(e/D) 4c 
c c-^ + (c + j 2 )(e/C) 4c ' 



= , / ,-2 



■Ji. 
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in agreement with Ref. JT^j. To lowest order in | — g, 
this one- loop result is consistent with Eq. (^|). In partic- 
ular, at the transition point j± = j z , the magnitude of 
the magnetization jump A = Ai(— 0) — M(+0) = 1 + j z 
agrees with our conjecture \f2g~ — (1 — j 2 ) _1 in the first 
order in j z . Approach to the zero-field value is logarith- 
mically slow at the transition point: A4(e) — A4(+0) oc 
1/lne for e > 0. It is worth mentioning that in the 
ferromagnetic case the magnetization is a nonuniversal 
quantity. In particular, models with different bandwidth 
cutoffs have values of M. different by <~ To be precise, 
our mapping is to the bosonized Kondo model (|Tl]). 

The magnetization jump in the ferromagnetic Kondo 
model is also related to the problem of a classical Ising 
chain with inverse-square interaction ]l6| ] , whose reduced 
Hamiltonian is (3 Hj = —K^^j SiSj/\i — j| 2 , Sj = ±1. 
The connection is established by expanding the partition 
function of the Kondo model in powers of Jj_, which is 
then identified with that of the Ising chain at K = l/Ag 
with additional short-ranged irrelevant interaction that 
depends on J± and bandwidth cutoff. The ID Ising 
model with 1/r 2 interaction has a finite-temperature 
phase transition fuifl . Below the transition temperature, 
or when K > K c , the Ising spins are ordered with a 
nonvanishing order parameter ^ = K(si} 2 . The order 
parameter vanishes discontinuously as the coupling K 
is reduced through the K c , a ID analog of the univer- 
sal jump at the Kosterlitz-Thoulcss transition. It has 
been shown rigorously that \E' > 1/2 in the ordered phase 
jl7j . Although unproved rigorously yet, it is believed that 
^ = 1/2 at K ~ K c , which has been confirmed recently 
by a large-scale numerical simulation [ p"8| |. This means 
that M(-0) - M(+0) = (sj) = \/2g, in agreement with 
our analysis. Note that the external parameter driving 
the phase transition is K in the Ising model whereas it 
is F/D in our model of a level coupled to a ID lead. 

So far we have ignored Coulomb repulsion between the 
resonant level and the ID lead. Assuming that the inter- 
action is short-ranged due to screening effects, we may 
model it by H u = u S z dip ^ ■ Since the unitary trans- 
formation gives U\fl u U2 = Hu + const, the effect of 
the Coulomb repulsion is to change the Kondo coupling, 
J z — > J z + 2iru. This gives rise to a shift of the phase 
boundary in Fig. [I] towards smaller g. In particular, at 
t — > the phase boundary is at g = ^(1 + t^)~ 2 - Physi- 
cally, this shift is due to Mahan's excitonic effect JlJ that 
enhances tunneling probability. 

It was recently demonstrated [[l9| |2(J that the charge 
of a quantum dot coupled to a Fermi-liquid lead can be 



measured with the aid of an electrometer based on a sin- 
gle electron transistor. To test our theoretical predictions 
one could perform a similar experiment with a lead in the 
Fractional Quantum Hall regime. A small dot coupled to 
the edge state in the lead would then play the role of the 
resonant level. In this setup one can vary the parameter 
g of the chiral Luttinger liquid in the edge state by ad- 
justing the magnetic field. The strength of coupling of 
the dot to the edge state in a GaAs heterostructure can 
be tuned by changing the appropriate gate voltage. 

In conclusion, we have shown that in the region of 
parameters shown in the phase diagram Fig. [I] the oc- 
cupation probability n(e) of a resonant level of energy 
e coupled to a Luttinger-liquid lead is a discontinuous 
function of e. The height of the discontinuity in n(e) at 
e = reaches ^/2g at the phase boundary. This picture is 
supported by mapping to a Kondo model, and interest- 
ing connections can also be made to the problems of the 
dissipative two-level system and the classical Ising chain 
with inverse-square interaction. 

We are grateful to S. Katsumoto for valuable discus- 
sions and to the Aspen Center for Physics where part of 
the work was performed. The work of AF was supported 
by the Grant-in- Aid for Scientific Research No. 11740199 
from the Japan Society for the Promotion of Science 
(JSPS). KAM acknowledges the support of the Sloan 
Foundation, the JSPS, and NSF Grant DMR-9974435. 



[1] For a review see, e.g., G. D. Mahan, Many-Particle 
Physics, 3rd ed. (Kluwer, New York, 2000). 

[2] A.M. Chang, L.N. Pfeiffer, and K.W. West, Phys. Rev. 
Lett. 77, 2538 (1996). 

[3] M. Bokrath et al, Nature (London) 397, 598 (1999). 

[4] Z. Yao et al, Nature (London) 402, 273 (1999). 

[5] O.M. Auslaender et al, Phys. Rev. Lett. 84, 1764 (1999). 

[6] M. Grayson et al, Phys. Rev. Lett. 86, 2645 (2001). 

[7] C.L. Kane and M.P.A. Fisher, Phys. Rev. B 46, 15233 

(1992) . 

[8] C. de C. Chamon and X.G. Wen, Phys. Rev. Lett. 70, 
2605 (1993). 

[9] A. Fumsaki and N. Nagaosa, Phys. Rev. B 47, 3827 

(1993) ; A. Furusaki, Phys. Rev. B 57, 7141 (1998). 
[10] X.-G. Wen, Int. J. Mod. Phys. B 6, 1711 (1992). 

[11] See, e.g., K.A. Matveev, Phys. Rev. B 51, 1743 (1995). 

[12] A.J. Leggett et al, Rev. Mod. Phys. 59, 1 (1987). 

[13] The second-order result ([]) is also found in P. Cedraschi 
and M. Biittiker, Ann. Phys. (NY.) 289, 1 (2001). 

[14] S. Chakravarty, Phys. Rev. Lett. 49, 681 (1982). 

[15] A.M. Tsvelick and PB. Wiegmann, Adv. Phys. 32, 453 
(1983); N. Andrei, K. Furuya, and J.H. Lowenstein, Rev. 
Mod. Phys. 55, 331 (1983). 

[16] PW. Anderson and G. Yuval, J. Phys. C 4, 607 (1971). 

[17] M. Aizenman, J.T. Chayes, L. Chayes, and CM. New- 
man, J. Stat. Phys. 50, 1 (1988); J.Z. Imbrie and CM. 
Newman, Commun. Math. Phys. 118, 303 (1988). 

[18] E. Luijten and H. Mefiingfeld, Phys. Rev. Lett. 86, 5305 
(2001). 



5 



[19] D. Berman, N.B. Zhitenev, R.C. Ashoori, and M. 

Shayegan, Phys. Rev. Lett. 82, 161 (1999). 
[20] D.S. Duncan et al, Appl. Phys. Lett. 74, 1045 (1999). 
[21] The Hamiltonian (^|) describing non-chiral bosons can 



be converted by a transformation <p(x) = y / 3^(l :E l) + 
sgn(x)6(\x\) / y/g to a chiral form similar to Eq. (H). 



